Newtonian Model Research Articles

THE ANALYSIS OF THE FLOW OF BLOOD IN A STENOSED ARTERY THROUGH SIMULATION: A COMPARISON AMONG VARIOUS NON-NEWTONIAN MODELS

This paper focuses on the dual quality of blood, Newtonian and non-Newtonian, in particular by exploring the energy curves. Careful investigation of the dual property of blood has been made by considering two different geometries to represent a stenosed arterial segment. We present a cautious assessment of non-Newtonian blood rheology impacts in arterial stream simulations by coupling the Newtonian and non-Newtonian models. The flow of energy through the two flow dimensions is meticulously investigated using velocity (kinetic energy), pressure, and wall shear stress (pressure energy). Besides, the proper implementation of an interface boundary condition (IBC) was emphasized to ensure consistency with the flow conditions downstream of a backward-facing step. The integration of the Newtonian and non-Newtonian models adjoins the novelty of the current research. The energy curves are obtained by implementing five different non-Newtonian models to designate a suitable non-Newtonian model for blood flow investigations. The combination of the non-Newtonian models enforced in this research is novel and particular attention is paid to the energy curves obtained. The conclusion was to elect the Carreau model as a suitable non-Newtonian rheological model for the blood flow study. This study was able to finalize the fact that the coupling of Newtonian and non-Newtonian models is necessary to obtain accurate results. For the sinusoidal waveform considered for the velocity, Carreau and the Power law models yield better results, eliminating the other non-Newtonian models from the list. With a better inlet condition imposed in the form of the Fourier series for pressure and velocity, the Carreau model yields the best results.

3D concrete printing using computational fluid dynamics: Modeling of material extrusion with slip boundaries

This paper investigates the role of slip boundary conditions in computational fluid dynamics modeling of material extrusion and layer deposition during 3D concrete printing. The mortar flow governed by the Navier-Stokes equations was simulated for two different slip boundary conditions at the extrusion nozzle wall: no-slip and free-slip. The simulations were conducted with two constitutive models: a generalized Newtonian fluid model and an elasto-viscoplastic fluid model. The cross-sectional shapes of up to three printed layers were compared to the experimental results from literature for different geometrical- and speed-ratios. The results reveal that employing free-slip boundary conditions at the extrusion nozzle wall improves layer-mimicking quality for both constitutive models, indicating the presence and importance of a lubricating layer of fine particles at the concrete-solid wall interface. This enhanced performance is primarily due to the observed decrease in extrusion pressure that minimizes layer height- and width-deviations compared to the experimental prints. Furthermore, the free-slip boundary conditions play an important role in predicting the multilayer prints, its deformation and groove shapes.

On the autonomous validation and comparison of particle models for a Newtonian laminar flow mixing model using PEPT

A methodology for the validation of a Computational Fluid Dynamics (CFD) mixing model using Positron Emission Particle Tracking (PEPT) was applied to a horizontal laminar stirred tank containing Newtonian glycerol. Particle trajectories recorded using PEPT were compared to three one-way particle-coupled CFD models (CFD-DEM, and CFD-Lagrangian, and CFD-Massless) through time-averaging and binning into a grid of 3-D voxels. A cell-by-cell comparison of particle velocities between the computational models and PEPT highlighted the accuracy of the CFD-DEM and CFD-Massless models over the Lagrangian model. Furthermore, using PEPT as a validation dataset the Pearson r-squared was used as a cost function to guide an autonomous calibration of a 5-dimensional DEM parameter space. A novel statistic, the simulation quality coefficient Qsim, was applied to assess accuracy between trajectory datasets and revealed the extent to which a model can be considered calibrated. Considering the balance between computational cost and precision, the CFD-massless model was the most effective particle model for this validation workflow. Further investigations will apply this methodology to calibrate non-Newtonian laminar CFD models.

Integrated drying model of lychee as a function of temperature and relative humidity

Drying is a universal method applied for food preservation. To date, several models have been developed to evaluate drying kinetics. In this study, lychee was dried employing a hot air dryer, and the drying kinetics was evaluated by comparing the Newtonian model, Henderson and Pabis model, Page model, and Logarithmic model. However, temperature and relative humidity, the key driving forces for drying kinetics, are not considered by these models. Thus, an integrated drying model, as a function of temperature and relative humidity, was developed to predict the hot air-drying kinetics and mass transfer phenomena of lychee followed by the calibration and validation of the model with independent experimental datasets. The model validation consisted of Nash- Sutcliffe model coefficient (E), coefficient of determination (R2) and index of agreement (d) and all of them were found close to 1 indicating perfect model fit. Besides, the developed model was applied for process optimization and scenario analysis. The drying rate constant was found as a function of temperature and relative humidity that was high at high temperature and low relative humidity. Interestingly, temperature showed a higher effect on the drying rate constant compared to relative humidity. Overall, the present study will open a new window to developing further drying model of lychee to optimize quality its quality parameters.

Applications of lubrication approximation theory in the analysis of the roll‐coating using a tangent hyperbolic fluid model

AbstractIn this study, the process of roll‐over‐web coating for a tangent hyperbolic fluid is theoretically examined. Fundamental laws of fluid mechanics are used to set up the flow equations and then simplified through lubrication approximation theory (LAT). The governing differential equations are solved with the help of the perturbation method in conjunction with the Runge–Kutta method. The outcomes of the study are presented through the graphs and tables. Study shows that enhancing the Weissenberg number upsurges the velocity and pressure while decreasing the pressure gradient. It is also noted that force shows an increasing trend while power is decreasing for bigger values . Results further depict that for a bigger value of the power law index , pressure gradient, pressure, force, and power decrease, and the coating thickness increases. It is worth mentioning here that for the shear thinning case, the tangent hyperbolic model predicts less pressure in the nip region when compared to the Newtonian model.

Stagnant points and uniform film analysis of Eyring fluid film flow on a vertically upward moving slippery flat plate

In this paper, we analyze how Eyring fluid film behavior on a vertically upward moving slippery flat plate can help predictive models in engineering, notably in coating and lubrication operations. This paper deals with the stagnant points and uniform film analysis of Eyring fluid film flow on a vertically upward moving slippery flat plate. The formulated ordinary differential equation is solved for exact analytic solution. Exact analytic expressions for velocity, flow rate, average velocity, shear stress components, and stagnant points are derived. A highly nonlinear algebraic equation is derived for film thickness. Equation is solved by Newton’s method using Maple code. The thickness of film widens with increasing relaxation time, constant plate velocity, and fluid density, while it decreases with increasing fluid viscosity and constant slip parameter. The analysis delineates that the positions of stagnant points tend to relocate closer the slippery plate as the Stokes number, Deborah number, and constant slip parameter increase. A high Deborah number enhances drainage, causing stagnant points to relocate closer to the slippery plate and the development of a stable elastic layer adjacent to the plate. For stagnant points and fluid film thickness, a comparison between the Eyring fluid and existing studies (EPTT, LPTT, UCM, and Newtonian fluid) is also provided. The validity of this paper with the existing studies is also presented by reducing Eyring fluid model to Newtonian model. The results of this research are significant for a wide range of biofluid applications, including agrochemical uses, paint and surface coating flow behavior, thin films on the cornea and lungs, and chemical and nuclear reactor design.

Accuracy and temporal analysis of non-Newtonian models of blood in the computational FFR – Numerical implementation

Non-invasive assessment of stenosed coronary is a goal for clinical applications. Thus, engineering aligned with medicine has been the solution to solve this problem. This work, serving as a proof of concept, has the main goal to predict the non-invasive fractional flow reserve (FFR) in an atherosclerotic coronary artery. A 3-element Windkessel and Womersley models were used to model boundary conditions of the patient-specific coronary arteries: the pressure at the outlets of the artery and the velocity at the inlet, respectively. Apart from implementing the previous conditions in the numerical code, the novelty of this work is to implement and use, in numerical simulations, Newtonian and non-Newtonian models for blood such as Carreau, Carreau-Yasuda, Casson and Simplified Phan-Thien/Tanner to analyse the accuracy of the computed FFR and the computational time of the simulations. Results indicate that the non-Newtonian viscoelastic model, Simplified Phan-Thien/Tanner, took the longest time but produced the most accurate results - 8h and 2.7% relative error between invasive and non-invasive FFR. The shear-thinning non-Newtonian models required lower computational times but returned a higher average error (3.9h; 7.6% error) which is considered significant. The Newtonian model required the least time but had the highest error value (2.9h; 7.8% error). Thus, although the computational time considering the non-Newtonian viscoelastic model is the highest, this model must be considered in the numerical simulations to obtain the FFR since it is the most accurate for this patient case. This research is a step forward to numerically quantify the non-invasive FFR of a patient.

Hemodynamic Study of Cerebral Arteriovenous Malformation: Newtonian and Non-Newtonian Blood Flow

ObjectiveArteriovenous malformation is a disease of the vascular system that occurs mainly in the cerebral arteries and spine. Numerical simulation as a powerful method is used to investigate the Cerebral Arteriovenous Malformation hemodynamic after occlusion of the abnormality step by step by embolization. MethodsThe CT Angiographic imaging data of two patients are used and a geometric model is extracted by the Mimics software. Numerical simulation of blood flow is performed in both Newtonian and non-Newtonian models. The Navier–Stokes, and continuity governing equations are solved by finite element method using the COMSOL Multiphysics software (the commercial computational fluid dynamics (CFD) simulation software). To validate the numerical results, the real data on blood flow rate in the feeding artery and draining veins are used, as well as angiographic images at different times. ResultsRegarding the comparison of pressure contours for different occlusions of 0, 30, 50, and 90%, by increasing the amount of occlusion in the nidus, there is an increase in the blood pressure. Regarding the comparison of the blood flow velocity in the feeding artery, draining veins, and inside the AVM nidus for Newtonian and non-Newtonian models, there is a significant difference between these two simulations in vessels with smaller dimensions (such as vessels inside the nidus). ConclusionBy increasing the amount of nidus occlusion, the blood pressure is increased, so the blood supply process is better. According to a significant difference between the Newtonian and non-Newtonian simulations in vessels with smaller dimensions (such as vessels inside the nidus), therefore, non-Newtonian simulation should be done for different occlusions of 30, 50, and 90%.

Oscillatory flow of rheological complex fluid in a flat channel

The article examines the oscillatory flow of rheological complex fluids in a flat channel, in which the rheological complex fluid was obtained based on the Maxwell model, and its mixture is described based on the Newtonian fluid model. Both mixtures of fluids are represented as a homogeneous model of a two-component liquid. In this case, differential equations of motion of homogeneous liquid mixtures are given. Based on this equation, the problem of oscillatory flows of rheological complex fluids in a flat channel is solved analytically. When solving, formulas are given for determining the longitudinal velocity. Using the obtained formulas, graphs of the longitudinal velocity distribution over the cross section of the channel are determined depending on the change in the oscillation frequency parameter, and with their help, appropriate conclusions are drawn.

Constraining the Milky Way Dynamic Mass and Newtonian Gravity through the Proper Motion of Messier’s Globular Clusters

The estimation of the proper motion of globular clusters is presented from the oldest observations of Messier's Globular clusters, comparing with NASA-NED and GAIA. The period-distance relationship, without assuming Dark matter models, allows us to limit the dynamic mass of the Milky Way and the invalidity of the Keplerian movement of Globular Clusters. The best fitting data are 14 2.3 0.3210 ± M for the dynamical mass, and β ≅ −5 3 in agreement with the generalization of the Virial Theorem in large-scale modified Newtonian gravity models.

Electroosmotic flow and heat transfer characteristics of a class of biofluids in microchannels at high Zeta potential

Peristalsis is an important dynamic phenomenon in the field of biomedical research, and has great application prospects in microscale fluids. In recent years, this biomimetic (peristaltic) phenomenon has gained widespread attention due to its large-scale applications in various medical and industrial fields, such as radiation therapy, peristaltic blood pumps, and drug delivery systems. In this study, the electroosmotic flow and heat transfer characteristics are investigated under high wall Zeta potential and slip boundary conditions for a certain type of biological fluid that satisfies the Newtonian fluid model. Fluid flows under the joint action of external electric field, magnetic field, and Joule heating. Firstly, without using the Debye-Hückel linear approximation, the numerical solutions are given by using the Chebyshev spectral method for the nonlinear Poisson-Boltzmann equation, the fourth-order differential equation satisfied by the stream function, and the thermal energy equation. The results are compared with those obtained by using the Debye-Hückel linear approximation to demonstrate the effectiveness of the numerical method used in this study. Secondly, the effects of wall Zeta potential, Hartmann number <inline-formula><tex-math id="M11">\begin{document}$H$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M11.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M11.png"/></alternatives></inline-formula>, electroosmotic parameter <inline-formula><tex-math id="M12">\begin{document}$m$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M12.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M12.png"/></alternatives></inline-formula>, slip parameter <inline-formula><tex-math id="M13">\begin{document}$\beta $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M13.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M13.png"/></alternatives></inline-formula> are discussed on the flow characteristics, peristaltic pumping, and trapping phenomena under electromagnetic environments, and the influence of Joule heating parameter <inline-formula><tex-math id="M14">\begin{document}$\gamma $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M14.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M14.png"/></alternatives></inline-formula> and Brinkman number <inline-formula><tex-math id="M15">\begin{document}$Br$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M15.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M15.png"/></alternatives></inline-formula> is explored on heat transfer characteristics. The results show that 1) wall Zeta potential plays an important role in controlling the velocity of fluid peristaltic flow; 2) the increase of electroosmotic parameter <inline-formula><tex-math id="M16">\begin{document}$m$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M16.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M16.png"/></alternatives></inline-formula> and slip parameter <inline-formula><tex-math id="M17">\begin{document}$\beta $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M17.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M17.png"/></alternatives></inline-formula> increases the flow velocity in the central region of the channel, while the increase of Hartmann number <inline-formula><tex-math id="M18">\begin{document}$H$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M18.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M18.png"/></alternatives></inline-formula> hinders the flow of fluid; 3) these flow behaviors exhibit opposite trends near the channel walls; 4) the number of streamlines captured by peristaltic transport decreases with Hartmann number <inline-formula><tex-math id="M19">\begin{document}$H$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M19.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M19.png"/></alternatives></inline-formula> and electroosmotic parameter <inline-formula><tex-math id="M20">\begin{document}$m$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M20.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M20.png"/></alternatives></inline-formula> increasing; 5) the increase of Joule heating parameter <inline-formula><tex-math id="M21">\begin{document}$\gamma $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M21.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M21.png"/></alternatives></inline-formula> and Brinkman number <inline-formula><tex-math id="M22">\begin{document}$Br$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M22.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20231685_M22.png"/></alternatives></inline-formula> leads temperature to rise.

Solving Numerically Ermakov-type Equation for Newtonian Cosmology Model with Vortex

Solving Numerically Ermakov-type Equation for Newtonian Cosmology Model with Vortex

Evolution of solitary hydroelastic strain waves in two coaxial cylindrical shells with the Schamel physical nonlinearity

The paper considers the formulation and solution of the hydroelasticity problem for studying wave processes in the system of two coaxial shells containing fluids in the annular gap between them and in the inner shell. We investigate the axisymmetric case for Kirchhoff–Lave type shells whose material obeys a physical law with a fractional exponent of the nonlinear term (Schamel nonlinearity). The dynamics of fluids in the shells is considered within the framework of the incompressible viscous Newtonian fluid model. The derivation of the Schamel nonlinear equations of shell dynamics makes it possible to develop a mathematical formulation of the problem, which includes the obtained equations, the dynamics equations of two shells, the fluid dynamics equations and the boundary conditions at the shell-fluid interfaces and at the flow symmetry axis. The asymptotic analysis of the problem is performed using perturbation techniques, and the system of two generalized Schamel equations is obtained. This system describes the evolution of nonlinear solitary hydroelastic strain waves in the coaxial shells filled with viscous fluids, taking into account the inertia of the fluid motion. In order to determine the fluid stress at the shell-fluid interfaces, we perform linearization of the fluid dynamics equations for fluids in the annular gap and in the inner shell. The linearized equations are solved by the iterative method. The inertial terms are excluded from the equations in the first iteration, while, in the second iteration, these are the values found in the first iteration. A numerical solution of the system of nonlinear evolution equations is obtained by applying a new difference scheme developed using the Gröbner basis technique. Computational experiments are performed to investigate the effect of fluid viscosity and the inertia of fluid motion in the shells on the wave process. In the absence of fluids in the inner shell, the results of calculations demonstrate that the strain waves in the shells during elastic interactions do not change their shape and amplitude, i.e., they are solitons. The presence of viscous fluid in the inner shell leads to attenuation of the wave process.

Cooling rate calculation of pure aluminum solidified during a novel gas-assisted continuous casting and extrusion process

Gas-assisted continuous casting and extrusion (GAC) process is a novel technique for the continuous preparation of non-equilibrium solidified materials. To evaluate the solidification rate of the melt during GAC, the cooling rate of pure aluminum melt solidified in the wheel groove is calculated, based on the Newtonian heat transfer model and the one-dimensional Fourier heat transfer model. Results show that the cooling rate of the melt during GAC mainly depends on the conduction heat transfer in the wheel groove and the melt thickness. When the thickness of the pure aluminum melt is between 1 and 2.5 mm, the cooling rate can reach 103–104 K s−1, which is much higher than that in CAXTEX (~10 K s−1). It means that GAC is a very effective process in the preparation of non-equilibrium solidification alloys.

Strong constraints on the gravitational law from Gaia DR3 wide binaries

ABSTRACT We test Milgromian dynamics (MOND) using wide binary stars (WBs) with separations of 2–30 kAU. Locally, the WB orbital velocity in MOND should exceed the Newtonian prediction by $\approx 20~{{\ \rm per\ cent}}$ at asymptotically large separations given the Galactic external field effect (EFE). We investigate this with a detailed statistical analysis of Gaia DR3 data on 8611 WBs within 250 pc of the Sun. Orbits are integrated in a rigorously calculated gravitational field that directly includes the EFE. We also allow line-of-sight contamination and undetected close binary companions to the stars in each WB. We interpolate between the Newtonian and Milgromian predictions using the parameter αgrav, with 0 indicating Newtonian gravity and 1 indicating MOND. Directly comparing the best Newtonian and Milgromian models reveals that Newtonian dynamics is preferred at 19σ confidence. Using a complementary Markov Chain Monte Carlo analysis, we find that $\alpha _{\rm {grav}} = -0.021^{+0.065}_{-0.045}$, which is fully consistent with Newtonian gravity but excludes MOND at 16σ confidence. This is in line with the similar result of Pittordis and Sutherland using a somewhat different sample selection and less thoroughly explored population model. We show that although our best-fitting model does not fully reproduce the observations, an overwhelmingly strong preference for Newtonian gravity remains in a considerable range of variations to our analysis. Adapting the MOND interpolating function to explain this result would cause tension with rotation curve constraints. We discuss the broader implications of our results in light of other works, concluding that MOND must be substantially modified on small scales to account for local WBs.

Predictions of the behavior of a single droplet and blends composed of Newtonian/viscoelastic minor phase and viscous major phase subjected to oscillatory shear flow

The oscillatory behavior of Newtonian and viscoelastic droplets in a Newtonian phase and blends composed of viscoelastic minor phase in a Newtonian major phase are theoretically investigated in this work. The non-Newtonian constrained volume model predictions are compared to experimental oscillatory shearing flow data of droplet and blends that are available in the literature. For a Newtonian droplet in a Newtonian phase, the model describes experimental droplet behavior well at high viscosity ratio and high strain amplitude. For viscoelastic droplet in a Newtonian phase, the model predicts less deformation for viscoelastic droplet than a comparable Newtonian droplet. For large amplitude oscillatory shear rheological data of blends composed of Boger fluid minor phase in a Newtonian major phase, the model shows improvement in prediction of the elastic modulus at high viscosity ratio, compared to the Newtonian model. The model also shows good agreement with large and minimum strain elastic moduli and large and minimum rate dynamic viscosities for small and large viscosity ratio viscoelastic polymer blends. For the blend of Boger fluid minor phase in a Newtonian major phase at viscosity ratio larger than one, we find that elasticity contributes to total stress from small to large strain amplitude values. For the blend of Boger fluid minor phase in a Newtonian major phase at viscosity ratio smaller than one, we find that elasticity is important only at large values of strain amplitude. Moreover, for the aforementioned blend at viscosity ratio larger than one, the predicted Lissajous Bowditch plots of excess stress do not reflect droplet shape/droplet orientation, and the opposite is true for the small viscosity ratio blend. Investigation of droplet long semi axis for the large viscosity ratio blend at LAOS conditions reveals oscillations at two to three times the imposed frequency.

Study on unsteady natural convection heat transfer of crude oil storage tank based on fractional-order Maxwell model

When considering the natural flow processes in crude oil in static storage tanks, researchers usually treat the nonlinear heat and mass transfer and flow using integer-order models. However, this ignores the nonlinear convection and heat-transfer phenomena that occur; these are dependent on time and space due to the viscoelasticity of crude oil, which is a non-Newtonian fluid. In this work, we examined the time and space dependence of the nonlinear natural convective heat-transfer processes in crude-oil storage tanks. Fractional-order Maxwell and Cattaneo–Vernotte models were applied to the traditional momentum and temperature equations considering the Boussinesq approximation. Based on the L1 discrete format, the fractional-order pressure-correction equation was derived, and an associated solver, CaputoTemporalFractionalFoam (CTfracFoam), was developed to simulate the convection processes in static crude-oil storage tanks. The results show that when α equals 0.1, the maximum velocity in monitoring area 2 is 0.1404 m/s, and the velocity gradually increases as α increases with a maximum value of 1.1468 m/s at 0.9. When the β is 0.1 and 0.2, the temperature changes in the transition region are −5.204 and −0.165 °C, respectively, while when β is greater than 0.2, the temperature changes are all positive values with a maximum value of 1.924 °C at 0.9. Compared with integer-order Fourier heat-transfer and Newtonian fluid-flow models, the developed fractional-order model can describe the natural convection phenomena of non-Newtonian fluids such as crude oil more flexible and accurately.

Prediction modeling of coastal sediment transport using accelerated smooth particle hydrodynamics approach

A GPU-accelerated 3D smooth particle hydrodynamics (SPH) scheme is developed and applied to a coastal multi-phase liquid-sediment interaction and sediment transport. The SPH scheme's meshless design and the sediment's particle structure enable the modeling of the waves' interactions with the sediment particles beyond the limitation of the mesh-based methods. A Newtonian constitutive model is used to model the liquid phase, and the sediment transport is formulated based on the Herschel-Bulkley-Papanastasiou (HBP) model. The yield characteristics of the sediment phase are estimated using the Drucker-Prager yield criterion. Due to the parallelization of the solution on graphics processing units, the 3D SPH scheme's performance, which uses millions of particles, is improved. Good correlations were observed in the SPH predictions and experimental measurements, with a maximum difference of 4.85 %. The validated scheme is applied to formulate forecasting models for the coastline sediment transport. It is found that erosion and scouring are expected at the coastline region inclined to the direction of the sea waves, with a predicted mass erosion of about 60e3 kg in four years. The wave's velocity is also established to be directly proportional to the sediment transport. The proposed multi-phase SPH methodology is proven effective for sediment transport prediction.

A Comparison of Newtonian and Non-Newtonian Models for Simulating Stenosis Development at the Bifurcation of the Carotid Artery

Blood is a shear-thinning non-Newtonian fluid in which the viscosity reduces with the shear rate. When simulating arterial flow, it is well established that the non-Newtonian nature is important in the smallest vessels; however, there is no consistent view as to whether it is required in larger arteries, such as the carotid. Here, we investigate the importance of incorporating a non-Newtonian model when applying a plaque deposition model which is based on near-wall local haemodynamic markers: the time-averaged near wall velocity and the ratio of the oscillatory shear index to the wall shear stress. In both cases the plaque deposition was similar between the Newtonian and non-Newtonian simulations, with the observed differences being no more significant than the differences between the selected markers. More significant differences were observed in the haemodynamic properties in the stenosed region, the most significant being that lower levels of near-wall reverse flow were observed for a non-Newtonian fluid.

Turbulent heat and mass lines in finite difference regime

AbstractFinite difference modelling of turbulent heat and mass flow visualization finds numerous applications in atmospheric flows/oceanic currents, wind turbines, thermal transfer in nuclear reactors, drag in oil pipelines, cooling of industrial machineries, and to investigate the complexity, dynamic and chaotic nature of the physical system. A turbulent phenomenon is effectively implemented in engineering, physics, earth sciences, bio‐engineering and medicine. Hence, motivated by the advantages of turbulence in various engineering fields, in the current article, a finite difference analysis is performed to demonstrate the k‐ε turbulence model‐based heat and mass lines visualization in boundary layer regime under turbulent buoyancy‐driven convective conditions along a cylinder. Turbulent flow characteristics are accurately explored by deploying the classical Newtonian flow model. Further, to accomplish a more sophisticated finite difference simulation, the effects of extra kinetic energy and its dissipation rate equations are considered. The produced Navier‐Stokes equations for time‐dependent turbulent heat and mass transmission are rendered to non‐dimensional by deploying suitable dimensionless numbers. The advanced coupled nonlinear turbulent unsteady buoyancy‐motivated vertical convection problem is then solved with a well‐sophisticated finite difference scheme such as Crank‐Nicolson technique using computational software. Authentication of current results with former solutions over a range of buoyancy number, Schmidt, and Prandtl parameters are presented. An extensive tabular and graphical discussion along with contours, heat and masslines visualization is included to enumerate the hydro‐dynamic, thermal and mass diffusion behaviour for the impact of emerged regulating numbers in the Prandtl regime. It is confirmed that, the accelerating turbulent buoyancy‐ratio number, maximizes the velocity, kinetic energy and dissipation rate at . Further, the numerical values of laminar thermal and mass diffusion rates are monotonically enhanced when compared to the turbulent values. Also, to verify the current findings, the authors compared the LRN k‐ε model turbulent results with the existing solutions and found good agreement. Further, the uniqueness and novelty of the current investigation is the exploration of heat and masslines in unsteady buoyancy‐driven convection regime under the influence of k‐ε turbulence model which extends the former studies and offers a more precise appraisal of the thermal and mass diffusion lines via the Crank‐Nicolson analysis.